Statistics of resonances and of delay times in quasiperiodic Schrodinger equations 
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We study the statistical distributions of the resonance widths P(r), and of delay times P(r) in 
one dimensional quasi-periodic tight-binding systems with one open channel. Both quantities are 
found to decay algebraically as F"", and t ' on small and large scales respectively. The exponents 
a, and 7 are related to the fractal dimension of the spectrum of the closed system as a = 1~\-D^ 
and 7 = 2 — Dq . Our results are verified for the Harper model at the metal-insulator transition and 
for Fibonacci lattices. 
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Quantum mechanical scattering, has been a subject of 
a rather intensive research activity during the last years. 
This interest was motivated by various areas of physics, 
ranging from nuclear [|][, atomic ||^ and molecular ||] 
physics, to mesoscopics^ and classical wave scattering 
1^. The most fundamental object characterising the pro- 
cess of quantum scattering is the unitary i'-matrix relat- 
ing the amplitudes of incoming waves to the amplitudes 
of outgoing waves. At present, there are two complemen- 
tary theoretical tools employed to calculate statistical 
properties of the 5— matrix, namely the semiclassical and 
the stochastic approach. The starting point of the first is 
a representation of the S*— matrix elements in terms of a 
sum over classical orbits while the later exploits the 
similarity with ensembles of Random Matrices (see 0] 
and references therein). Thus, for chaotic/ballistic sys- 
tems, many results are known. Among the most interest- 
ing is the knowledge of the Wigner delay time statistics 
and of the resonance width statistics . The former 
quantity captures the time-dependent aspects of quan- 
tum scattering. It can be interpreted as the typical 
time interval a scattering particle remains in the inter- 
action region. It is related to the energy derivative of 
the total phase shift ^{E) of the scattering matrix i.e. 
t{E) — The resonances represent long-lived inter- 

mediate states to which bound states of a closed system 
are converted due to coupling to continua. On a formal 
level, resonances show up as poles of the scattering ma- 
trix S{E) occurring at complex energies £„ = En — fTn, 
where E^ and r„ are called position and width of the 
resonances, respectively. 

Recently, the interest in quantum scattering has ex- 
tended to systems showing localization. For this case, 
there are analytical results about the distribution of 
phases of the S*— matrix and of delay times [p]-^. The 
former depends drastically on the disorder strength and 
energy |l2| , while for the latter a universal pow er law 
tail was found to hold |l^-|l^. Moreover, in a first 
analytical result about the distribution of resonances ap- 
peared. 

In the present paper we study delay time and reso- 



nance width statistics in a new setting, namely a class 
of systems, whose closed system analogues have fractal 
spectra. The latter exhibit energy level statistics that 
are in strong contrast to the level repulsion predicted by 
Random Matrix Theory (RMT) . Their level spacing 
distribution follows inverse power laws P{s) ~ which 
is a signature of level clustering. The power (3 was found 
to be related with the fractal dimension of the spectrum 
as (3 — 1 -\- Dq [|i^. Realizations of this class are, 
quasi-periodic systems with metal-insulator transition at 
some critical value of the on-site potential like the Harper 
model 1 16 1^, Fibonacci chains or quantum sys- 

tems with chaotic classical limit as the Kicked Harper 
Model jl9j. Here, for the first time we present conse- 
quences of the fractal nature of the spectrum in open 
systems. We consider open systems with one channel 
(the simplest possible scattering problem) and report the 
appearance of a new type of resonances width and delay 
time statistics. These distributions show inverse power 
law behaviour dictated by the fractal dimension Dq of 
the spectrum. Specifically, we show that the probabil- 
ity distributions of resonance widths 'P(r), and of delay 
times V{t) when generated over different energies, be- 
have as 
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For the calculation of 7'(r) and 7^(t) we employed two in- 
dependent approaches. Our results (|l|) are confirmed for 
two different types of quasi-periodic tight-binding models 
and are supported by analytical arguments. 

We consider a ID quasi-periodic sample of length L 
with one semi-infinite perfect lead attached on the left 
side. The system is described by the tight-binding Hamil- 
tonian: 

H^Y. + E d'^) + ii + I" + 1) ('^D (2) 

n n 

where Vn is the potential at site n. In the sequel we 
will consider examples where for < n < L, y„ is given 
by a quasi-periodic sequence. For n < 0, y„ = and 
we impose Dirichlet boundary conditions at the edge 
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tjjL+1 = 0. Therefore, for n < 0, scattering states of 
the form ■0n = e**"'" + 5'e~*'^" represent the superposi- 
tion of an incoming and a reflected plane wave. Here, 
k = arccos(£'/2) is the wave vector supported at the 
leads. Since there is only backscattering, the scattering 
matrix S{E) = e**^-^) is of unit modulus and the total in- 
formation about the scattering is contained in the phase 
One can write the scattering matrix in the form 

U2T 



S{E) EE e'*(^) = 1 - 2iw'^ sin k e 
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-e. (3) 



E — Hcs 

Ticff is an effective non-hermitian Hamiltonian given by 



(4) 



Hl is the part of the tight-binding Hamiltonian (g) with 
n = 0,....,L corresponding to the quasi-periodic sam- 
ple and e = (1, 0, 0, . . . , 0) ^ is an i— dimensional vector 
that describes at which site we couple the lead with our 
quasi-periodic sample. The strength of the coupling is 
given by w. In the sequel we will always consider w = I. 
Moreover, since arccos(i?/2) changes only slightly in the 
center of the band, we put E = and neglect the en- 
ergy dependence of TYcff- The poles of the S"— matrix are 
equal to the complex eigenvalues £ of Ticff- The latter are 
computed by direct diagonalization of Hcs- We note here 
that numerical diagonalization of complex non-hermitian 
matrices is a time consuming process and imposes limita- 
tions on the system size due to limited storage capacity. 
The size of the matrices that we used in our analysis 
below was up to rank 5000. 

For the calculation of the Wigner delay time r we have 
developed a simple iteration relation in 



in a unit cell of area a^. When a is an irrational number 
the period of the effective potential Vn is incommensurate 
with the lattice period. We consider generic irrationals 
which cannot be approximated "too well" by rationals. 
To this end we take a as the limit of successive rationals 
p/q, so that the potential becomes commensurate with 
the lattice with period q. Then we can define a scal- 
ing procedure where the incommensurate limit q —^ oo 
becomes equivalent with the thermodynamic limit. The 
states of the corresponding closed tight-binding system 
are extended when A < 2, and the spectrum consists 
of bands (ballistic regime). For A > 2 the spectrum is 
point-like and all states are exponentially localized (lo- 
calized regime). The most interesting case is the critical 
point A = 2 where we have a metal-insulator transition. 
At this point, the spectrum is a zero measure Cantor set 
with fractal dimension Dq < 0.5 j2^ while the states are 
critical, i.e. self-similar fluctuations of the wave function 
on all scales . 

First, we will investigate the statistical distribution of 
the resonance widths F, and delay times for the Harper 
model at the critical point A = 2. More exactly we de- 
termine the integrated distributions 
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sink J 1 + {tan{4>L - k) + AlY 
Gl = 1 + ALsin (2(0L - fc)) + Alcos\(j)L ~ k) 
tan((/)L+i) = tan{(l)L ~ k) + Al (5) 

where Al — VlI sinfc. Iteration relation (j^) has proved 
to be very convenient for numerical calculations since it 
anticipates the numerical differentiation which is a rather 
unstable operation. Moreover, it allows us to go to large 
system sizes. 

We motivate and numerically verify our results using 
first the well known Harper model which is a paradigm of 
quasi-periodic ID system with metal-insulator transition 

f,p7t. It is described by the tight-binding Hamiltonian 
with on-site potential given by 



Vn = Acos(27rcrn). 



(6) 



This system effectively describes a particle in a two- 
dimensional periodic potential in a uniform magnetic 
field with a = a^eB/hc being the number of flux quanta 



Vint (x) 



V{x)dx 



(7) 



whose derivatives V{x) = —dVint/dx determine the 
probability density of resonance widths V{x = F) and 
delay times 'P{x ^ t). In all our calculations we will 
take approximants of the golden mean ctg = (\/5 — 
For this case it is known that Dq w 0.5 [ pOj . 

Figure 1 shows T'mt(r) for two different rational ap- 
proximants a of the golden mean aa. It clearly displays 
an inverse power law 



^™t(r) ^Fi 



(8) 



and thus the resonance width distribution behaves as 
stated in (|^) with a ~ 1.5 = 1 -I- Dq . The integrated 
resonance width distribution cuts off at a small value of 
F's (see Fig. 1), since for all rational approximants of 
ac the total number of is finite. This cutoff, how- 
ever, can be shifted to arbitrarily small values for higher 
approximants. 

Next we investigated the delay time statistics 7'(t). In 
Fig. 2 we report the integrated Vint{T) for three differ- 
ent rational approximants of the golden mean. Due to 
the efficiency of our iteration relation (^ we can approx- 
imate <7g by increasing the periodicity q of the potential 
as much as we like. Our numerical data are in agreement 
with an inverse power law i.e. 



V^ntiT) 



(9) 



with a value of7~1.5 = 2 — given by a best least 
square fit, in perfect agreement with Eqn. (m). 
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The connection between the exponents a, 7 and the 
fractal dimension of the close system calls for an 
argument for its explanation. The following heuristic ar- 
gument, similar in spirit to [ p^ , p2[ provides some under- 
standing of the power laws (□). We consider successive 
rational approximants Ui ~ Villi of the continued frac- 
tion expansion of a. On a length scale qi the periodicity 
of the potential is not manifest and we may assume that 
a wave packet moves as var{t) ~ t'^^o We attach 

the lead at the end of the segment qt which results in 
broadening the energy levels by a width T. The max- 
imum time needed for a particle to recognize the exis- 

fence of the leads, is Tq. ^ q^ " . The latter is related 
to the minimum level width Tq. ^ 1 /xq. . The number of 
states living in the interval is ~ qi and thus determines 
the number of states with resonance widths T > l/r^^. 
Thus Vint{rqi) ~ qi ^ r^-^o , By repeating the same 
argument for higher approximants di+i = pi+i/g^+i we 
conclude that V{T) ~ r^(^+^o \ in agreement with (|^). 
Although the numerical results support the validity of 
the above argument, a rigorous mathematical proof is 
still lacking. 

Next, we present another argument, which allows us to 
understand the relation between the power law decay ex- 
ponent 7 and the fractal dimension Dq i.e. j — 2 - 
Our starting point is the well known relation 



;^^(i?-£;„)2 + r2/4 



(10) 



which connects the Wigner delay times and the poles 
of the 5— matrix. It is evident that anomalously large 
time delay t{E) ^ F^^ corresponds to the cases when 
E ~ En and r„ ^ 1. In the neighbourhood of these 
points, t{E) can be approximated by a single Lorentzian 
(|lO|). Samphng the energies E with step AE <C Tmin 
we calculate the number of points for which the time 
delay is larger than some fixed value r. Assuming that 
the contribution of each Lorentzian is proportional to its 
width one can estimate this number as ^i/^. Tn/ A.E. 
For the integrated distribution of delay times we obtain 
V^nt{T) ~ f^"' drp(T)T ~ r-(2-") in the limit A£; ^ 
where we used the small resonance width asymptotics 
iyen by Eqn. (|l]) (for similar argumentation see also 
5|). Then for the asymptotic distribution of delay 
times we get V{t) ~ r~(^~^o ^ in agreement with (|^) 
and our numerical findings. 

The validity of the heuristic arguments (and thus of 
Eqs. (|^)) can be verified in more cases in the Fibonacci 
chain model of a one dimensional quasi-crystal where 
other scaling exponents can be obtained. Here the poten- 
tial Vn only takes the two values +V and —V arranged 
in a Fibonacci sequence . It was shown that the spec- 
trum is a Cantor set with zero Lebesgue- measure for all 
y > 0. We again find inverse power laws for the inte- 



grated distributions 'P(r) and ■P(r). Here the exponent 
depends on the potential strength V, while Eqs. (P still 
relate the corresponding statistics to the fractal dimen- 
sion Dq . Our results for various V values are summa- 
rized in Fig. 3 and show a nice agreement between the 
exponents a, 7 and according to Eq. (|]). 

Because of lack of space we defer the discussion of other 
results, like the fractal nature of the resonance widths, 
and the behaviour of the delay time autocorrelation func- 
tion to a later publication |^^ . 
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FIG. 1. Vinti^) of the Harper model (A = 2) for three approx- 
imants of ctg , cti = .j^: ; (72 = ^||| ; and 0-3 = |y|y ■ An inverse 
power law Pint{r) ^ t^~" is evident. A least squares fit yields 
a pa 1.5 in accordance with Dq ~ 0.5 and Eqn. (|l|). As is seen the 
lower cutoff of the scaling region decreases for higher approximants. 




FIG. 2. VintiT) of the Harper model (A = 2) for three approxi 
mants of the golden mean cti = ||| ; 02 = ; and 03 — ^^^O'"' 
An inverse power law Pint{T) ~ '' is evident 
fit yields 7 ^ 1.5 in accordance with ~ 0.5 and Eqn. ( 
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is seen the upper cutoff of the scaling region increases for higher 
approximants . 
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FIG. 3. Power law exponents a, 7 (plotted as a — 1 and 2 — 7) of 
the resonance widths and of the delay time distributions, respec- 
tively, as a function of the potential strength V for the Fibonacci 
model. We also plot the fractal dimension Dq of the spectrum 
(the solid line is to guide the eye). Our numerical data show that 
a and 7 are related to the Hausdorff dimension Dq" according to 
Eqns. (Fl). 
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